ON THE ARTIFICIAL COMPRESSIBILITY METHOD FOR 
THE NAVIER STOKES FOURIER SYSTEM 

DONATELLA DONATELLI 

Abstract. This paper deals with the approximation of the weak solu- 
tions of the incompressible Navier Stokes Fourier system. In particular 
it extends the artificial compressibility method for the Leray weak so- 
lutions of the Navier Stokes equation, used by Temam [31], in the case 
of a bounded domain and later in [5] in the case of the whole space. 
By exploiting the wave equation structure of the pressure of the ap- 
proximating system the convergence of the approximating sequences is 
achieved by means of dispersive estimate of Strichartz type. It will be 
_^ ■ proved that the projection of the approximating velocity fields on the 

fvj ' divergence free vectors is relatively compact and converges to a weak 

solution of the incompressible Navier Stokes Fourier system. 
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This paper is concerned with the artificial compressibiHty approximation 
for the following Navier Stokes Fourier system, namely 

' dtu + {u • \7)u — fiAu = Vp 



(1.1) 
div n = 

u{x,Q) = uq{x), 6'(x,0) = 6*0(3;) 
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[~>^ . where (x, t) e M^ x [0, T], u G M^ denotes the velocity vector field , p G M the 

CD \ pressure of the fluid , E M is the temperature of the fluid, ^ is the kinematic 

viscosity and k is the heat conductivity. The model (|l.ip was proposed in 

[12j and [TU] and, as we can observe the system (jl.ip is an extension of 

the Navier Stokes equation which governs both the velocity fluid u and 

S^ . the fluctuation of the temperature fleld 6 in an incompressible fluid. Let 

H \ us point out that in the model under consideration the temperature fleld 

is just advected by the velocity fleld u and diffuses according to Fourier's 

law. More complicated effects such as viscous heating for example (see |19j ) 

will be not considered in the present paper, but they will be treated in 

future work. The artiflcial compressibility approximation was introduced 

by Chorin [U [2] , Temam [29l |30] and Oskolkov [21] , in order to deal with 

the difflculty induced by the incompressibility constraints in the numerical 

approximations to the Navier Stokes equation. The paper of Temam [291 

[30] and his book [3l] discuss the convergence of these approximations on 
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bounded domains by using the classical Sobolev compactness embedding 
theorems and they recover the compactness in time by the classical Lions 
|18j method of fractional derivatives. In [5] the authors cannot make use 
of the classical compactness theorems, therefore they exploit the underlying 
wave equation structure and recover the necessary compactness by means of 
dispersive type estimates. Here it will be used a similar approach adapted 
to the artificial compressibility approximation of the Navier Stokes Fourier 
system. In particular the following system we will considered 



(1.2) 



where (x,t) G M^ x [0,T], n^ = u^ix^t) G M^ r = r(x,t) G M and 
p^ =ps(x,f) G M. 

The system will be regarded as a semilinear wave type equation for the 
pressure function and the dispersive estimates will be carried out by using 
certain classical L^-type estimates due to Strichartz [HI [I5l [28] . The partic- 
ular type of Strichartz estimates that we are going to use here can be found 
in the book of Sogge [26] or deduced by the so called bilinear estimates of 
Klainerman and Machedon [16j and Foschi Klainerman [Oj. Our analysis can 
also be related to the convergence of the incompressible limit problem via 
a formal expansion (see for instance Temam [ST, Chapter 3). In particular 
our wave equation structure has some similarities with the one exploited in 
various way by the paper of Desjardin, Grenier, Lions, Masmoudi [1], Des- 
jardin and Grenier [3]. 

The interest into the artificial compressibility methods started with the pre- 
viously mentioned results of Chorin and Temam and was later on inves- 
tigated by Ghidaglia and Temam [10|. Later developments of numerical 
investigations in the directions of projections methods have been carried 
out by [H, [8], [22], [23], [20], [T7], [25j, [32j. 

The approximation of the Navier Stokes Fourier system has also been ana- 
lyzed by a different point of view. In fact F. Golse and L. Saint-Raimond 
in [12] , [13] recover weak solutions for the system (jl.ip in the frame of the 
hydrodynamic limit for the Boltzmann equation. 

This paper is organized as follows. In Section 2 we recall the mathematical 
tools needed in the paper and we recall same basic definitions. In Section 3 
we set up our problem, we explain our approximating system and we state 
our main result. The Section 4 is devoted to recover the a priori estimates 
needed to get the strong convergence of the approximating sequences and 
to prove the main theorem. Finally, in Section 5, we give the proof of the 
main result. 

2. Preliminaries 

For convenience of the reader we establish some notations and recall some 
basic facts that will be useful in the sequel. 

We will denote by V{W^ x R+) the space of test function C^{W^ x M+), 
by P'(M'^ X M+) the space of Schwartz distributions and (•,•) the duality 
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bracket between V and V and by MtX' the space C°([0, T]]X)'. Moreover 
W^'P{W^) = {I - A)-^LP{R'^) and H''{R'^) = W'''^{R'^) denote the nonho- 
mogeneous Sobolev spaces for any 1 < p < oo and fc G M. W '^{M. ) = 
{-Ay^LP{M.'^) and H''{M.'^) = Ty'='2(M'^) denote the homogeneous Sobolev 
spaces. The notations L^L% and L^Wx'"^ will abbreviate respectively the 
spaces LP([0,r];L'?(M'^)), and LP{[0,T];W''''i{R'^)). 

The operators Q and P denote the Leray's projectors on the space of gradi- 
ents vector fields and on the space of divergence - free vector fields, respec- 
tively. Hence, in the sense of distribution, one has 

Q = VA"Miv P = I-Q. (2.1) 

Let us remark that Q and P can be expressed in terms of Riesz multipliers, 
therefore they are bounded linear operators on every W 'P (1 < p < cxd) 
space (see 



Let us recall that if w is a (weak) solution of the following wave equation in 
the space [0, T] x R'^ 

^ + A)w{t,x)=F{t,x) 



for some data f,g,F and time < T < oo, then w satifies the following 
Strichartz estimates, (see [H], [T5] ) 

lkllL?Lj + \\dtw\\^,^-i,r < 11/11^. + ||5||^.-i + \\F\\^,'^,,, (2.2) 

where {q,r), {q,r) are wave admissible pairs, namely they satisfy 
2 . , ^ /I 1\ 2 , , , /I 1 

and moreover the following conditions holds 

1 d d Id 

- + - = 77-7 = - + -- 2. 
q r 2 q' r 

Later on we shall use (12. 2|) in the case of d = 3, {q',r') = (1,3/2), so one 
has 7 = 1/2 and {q,r) = (4,4), namely the following estimate holds 

II^IIl^, + \\dtw\\^,^-iA < 11/11^1/2 + ll5ll^-i/2 + 11^11^14/2- (2-3) 

Beside the Strichartz estimate (12.21) or (12. 3|) in the case of d = 3 (see [26]). 
there exists a non standard estimate which follows from an earlier linear 
Strichartz ^E\ estimate. This inequality can also be deduced by the bilinear 
estimates of Klainerman and Machedon [16], Foschi and Klainerman [2], 
namely 

\MlI^ + \\dtw\\^,^-iA < 11/11^1/2 + \\g\\^i/2 + WFhiLl- (2.4) 

Finally, we mention here the following lemma that will be usefull in the next 
sections, 
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Lemma 2.1. Let us consider a smoothing kernel j G Cq°(M'^), such that 
j > 0, Jjgd jdx = 1, and define the Friedrichs moHifiers as 

ja{x) = a-'^j (^^) . 

Then for any f G H^{M. ), one has 

ll/-/*Ja|lLp(R^) <Cpa^''^||V/||^2(K.), (2.5) 

where 

p£[2,oo) ifd = 2, p£[2,Q] if d = 3 and a = d( 

V2 p 

Moreover the following Young type inequality holds 

\\f*Ja\\Lp(m < Ca''''^^''^h\f\\w-'^'>(m, (2-6) 



for any p,q £ [1, oo], q < p, s > 0, a £ {0, 1). 

3. Approximating system and main result 

Let us consider the incompressible Navier Stokes Fourier system (jl.ip we 
recall (see P. L. Lions |19j ) the notion of weak solution that will be used later 
for (fTTTl . 

Definition 3.1. We say that the couple (n, 6) is a weak solution of the Navier 
Stokes Fourier system if it satisfies 

rT r f f)^n\ r 

uq ■ ipdx, 



n( uVn • V(i? — UiUjdiifj — u ■ -—- ] dxdt = j 
_ _; A dt J J^c 



kV9 ■ Vip — On ■ Vip — 9 ■ -^ ] dxdt = / Oq ■ tpdx, 
/O JR'i V ot 




and 



for all if G C^{R'^ x [0,T]), div^j = and 

div'u = in2?'(IR'^x [0,T]) 
Moreover the following energy inequality holds 

1 r '•* 



{\u{x,t)\^ + \e{x,t)\^)dx+ / / {fi\Vu{x,s)\' + K\V9{x,s)\^)dxds 

<o / {\uo\^ + \Oof)dx, forant>0. (3.1) 

Sometimes one refers also to this solution as Leray weak solutions al- 
though Leray himself didn't study thermal effects. It is worth to mention 
that the conservation of energy (j3.1|) is strictly related to the space dimen- 
sion d = 3 because of the integrability requirements of |np or u9. For a 
detailed mathematical description of models involving temperature and on 
the connected mathematical difficulties in studying them we refer to the 
book of P.L. Lions [19]. In order to approximate the svstem (jl.ip we wish to 
use the system (jl.2p . As we can observe the first two equations of the system 
(jl.2p resemble respectively the equation for conservation of momentum and 
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the equation of the temperature balance for the system (jl.ip except for the 
two terms 

We need to add these two correction in order to avoid the paradox of in- 
creasing the kinetic energy along the motion. In the third equation of the 
system ()1.2p we introduce a "linearized" compressibility constraint given by 
the equation 

dtp^ = — div u^ . 

e 

The limiting behaviour as e J, of the initial data of (|1.2|) deserves a little 
discussion. Indeed to solve p.2|) . three initial conditions are required 

u'{x,^)=ul{x), r(x,O)=0g(x), p'{x,Q)=pl{x), (3.2) 

while the Navier Stokes Fourier system (jl.ip requires only an initial condi- 
tion on the velocity u and an initial condition on the temperature 9. Hence 
our approximation will be consistent if the initial datum on the pressure 
will be eliminated by an "initial layer" phenomenon. Since in the limit we 
have to deal with weak solutions that verify (|3.ip it is reasonable to require 
the finite energy constraint to be satisfied by the approximating sequences 
{u^,9^,p'^). So we can deduce a natural behaviour to be imposed on the 
initial data {uq,9q,Pq), namely 

ug = u%-,0) — >uo = n(-, 0) strongly in L^{M.^) (ID) 

91, = u%-,0) — >9o = 9{-,0) strongly in L'^{R^) 

^epl = Vep^-, 0) — >0 strongly in L^{W^). 

Let us remark that the convergence of y/epQ to is necessary to avoid the 
presence of concentrations of energy in the limit. 

Since it will not affect our approximation process, for semplicity from now 
on, we will take n = 1 and k = 1. We will omit here the proof of existence 
of solutions for the system ()1.2|) since it can be obtained via standard fi- 
nite dimensional Galerkin type approximations. For the convenience of the 
reader we summarize here this existence theorem. 

Theorem 3.2. Let {uq,9q,Pq) satisfy the conditions (ID), for some e > 0. 
Then the system (|1.2p has a weak solution {u'^,9^,p^), with the following 
properties 

(i) u" GL°°([0,r];L2(R3))nL2([0,r];iJi( 

(ii) re L°°([o,r];L2(M3))nL2([o,r];ii'M 

(ill) ^p'eL^i[0,T]-L^iR^)), 
for all T > 0. 

Let us now state our main result. The convergence of {u^} will be de- 
scribed by analyzing the convergence of the associated Hodge decomposition. 

Theorem 3.3. Let {u^,9^,p'') be a sequence of weak solution in W^ of the 
system ()1.2p . assume that the initial data satisfy (ID). Then 

(i) There exists u £ L'=^{[0,T]; L^{R^)) n L^{[0,T]; H^{R^)) such that 
u^^u weakly in L'^{[0,T];II^{R^)). 
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(ii) The gradient component Qu^ of the vector field u^ satisfies 

Qu^ — > strongly in L^([0,r]; ^^(R^)), for any p e [4,6). 

(iii) The divergence free component Pu'^ of the vector field u^ satisfies 

Pu" — >Pu = u strongly in L^{[Q,T]]LJ^^{^^)). 

(iv) There exists 9 G L°°([0,T]; L2(IR3)) p L2([0,r]; 7^1(^3)) such that 

r — >e strongly in L^{[{), T];Lf^^{m?)), 

Vr ^ Ve weakly in L^{[d, T] x R^). 

(v) The sequence {p^} will converge in the sense of distribution to 

p = A-^ div {{u ■ V)u) = A'Hr{{Duf). 

(vi) u = Pu and 9 are weak solutions to the incompressible Navier Stokes 
Fourier system in the sense of the Definition (j3.ip and the following 
energy inequality holds for all t G [0, T], 



{\u{x,t)\'^ + \e{x,t)\'^)dx+ I I {\Vu{x,s)\^ + \Ve{x,s)\^)dxds 

i JO JRS 



<2 1 i\Mxr + \ooixr)dx. 



4. A PRIORI ESTIMATES 



In this section we wish to estabhsh the a priori estimates, independent 
on £, for the solutions of the system (II. 2p which are necessary to prove the 
Theorem 13.31 We will achieve this goal in two steps. First of all we will 
recover the a priori estimates that come from the classical energy estimates 
related to the system p.2p . Then, we get stronger estimates by exploiting 
the structure of the system. In fact, as we will see later on, the sequence p^ 
satifies a wave type equation. This will allow us to apply to p'^ the Strichartz 
estimates (j2.2p . (|2.4p and to get in this way dispersive bounds on p^. 



4.1. Energy estimates. The next results concerns the energy type esti- 
mate for the system (|1.2p . 

Theorem 4.1. Let us consider the solution {u^,0^,p^) of the Cauchy prob- 
lem for the system (jl.2p . Assume that the hypotheses (ID) hold, then one 
has 

E{t) +11 (|Vn^(x, s)|2 + |vr (x, s)\'^)dxds = E{0), (4.1) 

where we set 

Eit) = / (^|u^(x,t)p + ^|r(x,t)p + '-\f{x,t)A dx. (4.2) 



Proof. We multiply, as usual, the first and second equation of the system 
()1.2p respectively by u^ and 0^ and the third by p^ , then we sum up and 
integrate by parts in space and time, hence we get (14. ip . El 
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Corollary 4.2. Let us consider the solution {u'^,6^,p^) of the Cauchy prob- 
lem for the system (jl.2p . Let us assume that the hypotheses (ID) hold, then 
it follows 

yfep^ is bounded in L°°{[0,T];L'^{R^)), (4.3) 

epf is relatively compact in H^'^{[0,T] x R^), (4.4) 

VM^ Vr are bounded in L^{[0, T] x R^), (4.5) 

M^r are bounded in L'^{[0,T]; L'^{R^)) n L'^{[0,T];L^{R^)), 

(4.6) 



^^^eylf are bounded in L'^{[0,T];L^{M.^))nL^{[0,T];L^/^{R^)), 



(4.7) 

(div?)r "^' ^""^^'^ ^" L2([0,r];Li(R3))nLi([0,r];L3/2(R3)). 

(4.8) 

Proo/. (113]), (fO]) . (|i3]) follow from (fil^ . while ([32D follows from (|iT^ and 
Sobolev's embeddings theorems. Finally (j4.7p and (j4.8p come from (j4.5p . 
CED. □ 



4.2. Pressure wave equation. In this section we wish to get stronger 
bounds on p^. First of all, we observe that p^ satisfies a wave equation. In 
fact, by differentiating with respect to time the equation (jl.2P o and by using 



(|1.2p ^ we have 

edttp^ - Ap^ + A div u^ - div ( (u^ • V)u^ + - (div u^)u^] =0. 

Now, by rescaling the time variable as 

t 

and consequently the velocity, the temperature and the pressure in the fol- 
lowing way 

u{x,t) = u^{x,y/eT), 9{x,t) =9'^{x,y/eT) p{x,t) =p'^{x,y/eT), (4.9) 

we get that p satisfies the following wave equation 

5rrP- Ap + A div u- div ((n-V)n + -(div n)ttj =0. (4.10) 

This structure allow us to use on p the Strichartz estimates (|2.4p , ()2.2p . Now 
the analysis for p will follow the same line of arguments as in [5] . So we split 
p as p = pi + p2 where pi and p2 solve the following wave equations: 

fSrrPi - Api = -Adivu = Fi 
ypi {x, 0) = drPi (x, 0) = 0, 

drrP2 - Ap2 = div f (u • V) u + -(div u)u j = F2 

^P2{X,0) =p{x,0) drP2{x,0) =drP{x,0). 

Therefore we are able to prove the following theorem. 
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Theorem 4.3. Let us consider the solution {u^,0^,p^) of the Cauchy prob- 
lem for the system (j2.ip . Assume that the hypotheses (ID) hold. Then we 
set the following estimate 

+ VTll divn'^ll r2r2 



1 

(4.13) 



+ II {u" • V)u^ + ^(divii^)ii^ii ^1^3/2 



Proof. Since pi and p2 are solutions of the wave equations (I4.11|) . (I4.12|) . 
we can apply the Strichartz estimates ()2.3p and (12. 4p . with (x,r) G M^ x 
{0,T/^/e). Since A~^pi satisfies the equation 

drr{A~'pi) - A(A-ipi) = A~iFi, (4.14) 

then by using the Strichartz estimates ()2.4p we get 



ll^illi^iy-^-* + \\9rPi\\L4w-^A < ^11 dWuW^L^. (4.15) 

In the same way we have that A~^''^p2 satisfies the equation 

drAA-'/^P2) - A(A-V2p,) = A-^/^F2, (4.16) 

therefore by using the estimate p.Sp we obtain 

||P2|li4iy-i.4 + \\drP2\\Liw-^'^ ~ IIpC^:, 0) || -1/2 + ||9^p(a;, 0) || 3/2 

+ \\{u-V)u+-{divu)u\\.^.3/2 (4.17) 

Now by using (j4.15p . (j4.17p it follows that p verifies 

MlfW-'^^ + W^tPWl^^w-''^ ^ II:PiIIl4H/-2-* + IIP2|li4H/-M (4.18) 

+ \\9tPi\\l4w-3A + \\drP2\\L4w-2A 
< \\p{x,0)\\^-l/2 + \\drP{x,0)\\^-S/2 

Vt 1 

+ -JJI II *^^^ ^Wl^lI, + II (n • V) n + -(div u)u\\^^^^. 
Finally, since 

II^IIl'-((0,T/v/^);L9{M3)) = £ ''IIP^IIl'-([0,T];L'J(]R3)) 

we end up with ()4.13p . D 

5. Strong convergence 

In this section we conclude the proof of the Theorem 13.31 We will show 
the strong convergence of Qu^, Pu^ and 6^. In particular we will be able 
to prove that the gradient part of the velocity Qu'^ converges strongly to 
0, while the incompressible component of the velocity field Pu^ converges 
strongly to Pu = u, where u is the limit profile as e J, of n^. Moreover 
we will be able to prove some time regularity of 9^ that entails the strong 
compactness for 9^ and consequently the strong convergence to 9, solution 
of the system (jl.ip . 
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5.1. Strong convergence of Qu"^, Pu'^ and 6^. We start this section with 
some easy consequences of the a priori estimates estabhshed in the previous 
section. 

Proposition 5.1. Let us consider the solution {u^,6'^,p^) of the Cauchy 
problem for the system (jl.2p . Assume that the hypotheses (ID) hold. Then, 
as e I 0, one has 

ep' — ^0 strongly in L°°([0,r]; ^^(M^)) p L^([0,T]; ^-^'^(M^)), 

(5.1) 

div-u" — >0 siron^/y in M/^-^'°°([0,r];L2(IR3))nL^([0,r];VF-3'^(M^)). 

(5.2) 

Proof. dSl]), (1521) follow from the estimates ([O]), (liTl and the second 
equation of the system ()2.ip . El 



Now, we wish to show that the gradient part of the velocity field Qu^ goes 
strongly to as e | 0. As we will see in the next proposition, this will be a 
consequence of the estimate (|4.13|) and of the Lemma 12.11 

Proposition 5.2. Let us consider the solution {u^,9'^,p^) of the Cauchy 
problem for the system, (|1.2|) . Assum,e that the hypotheses (ID) hold. Then 
as e [Q, 

Qu' — > strongly in L'^{[0, T];LP{R^)) for any p £ [4, 6) . (5.3) 



Proof. In order to prove the Proposition 15.21 we split Qu^ as follows 

where ja is the smoothing kernel defined in Lemma 12.11 Now we estimate 
separately Ji and J2. For Ji by using ()2.5p we get 



1-3(1-^)^ r„V7^.,.,.M,2 ^.A'^' l-3(^ ^^ 



Ji<a "V2 v) \ \\VQu^{t)\\i2dt\ <a '\^ ''^ llVn^lLa^a. (5.4) 

Hence from the identity Qu"^ = —e^'^'S/A~^e'^'^dtp and by the inequality 
(|2.6p we get J2 satisfies the following estimate 



J2 < e'/'a-'~'^-^~-^)T'/%'/%p\\^,^-^,,. (5.5) 

Therefore, summing up (j5.4p and ()5.5p and by using ()4.5p and (|4.13p . we 
conclude for any p £ [4, 6) that 



1 M -, ,„ o 0/ 1 1 ' 



\\Qu'\\L2rP < Ca'-'y^'^) +CTe''^a ^'■'K^'pJ . (5.6) 

t ^ 

Finally we choose a in terms of e in order that the two terms in the right 
hand side of the previous inequality have the same order, namely 

a = e^/i^ (5.7) 

Therefore we obtain 



WQu'Wl^lp < Cts'^^p for any p G [4, 6). 



D 



'{t + h)-e^t)\\l2,,nT]xm= / dtdxw' ■ {w' - w' * jc,) 
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It remains to prove the strong compactness of the incompressible compo- 
nent of the velocity field and of the temperature. To achieve this goal we 
need to recall here, the following theorem (see |24] ) . 

Theorem 5.3. Let be J^ C LP{[0,T];B), 1 <p <oo, B a Banach space. T 
is relatively compact in LP{[0,T]; B) for 1 < p < oo, or in C{[0,T]; B) for 
p = oo if and only if 

(i) < / f{t)dt, / G i? > is relatively compact in B , < ti < t2 < T, 

(ii) lim \\f{t + h) - f{t)\\LP{[o,T-h];B) = uniformly for any f &T. 

The compactness can be obtained by looking at some time regularity 
properties of Pu^ and Q^ and by using the Theorem 15. 3t but before we need 
to prove the following lemmas. 

Lemma 5.4. Let us consider the solution {u^,9^,p^) of the Cauchy problem 
for the system (jl.2p . Assume that the hypotheses (ID) hold. Then for all 
h G (0, 1), we have 

||r(t + /i) -r(t)||i2([o,T]xR3) < Crh'/'. (5.8) 

Proof. Let us set w^ = 9^{t + h) — 9^{t), we have 

2([0,T]xM3) - J^ 

+ dtdxw" ■ {w^ * ja) =Il + l2- (5.9) 

Jo Jm.3 

By using (j2.5p we can estimate /i in the following way 

Ii < Ww'Wl^lI / \\w'{t) - {w' * ja){t)\\L2dt 

<ari/2|m£|| ||Y^£|| _ (5_;lo) 

r^ II lli^t i^^ll 117.^,^ V J 

Let us reformulate w'^ in integral form by using the equation (ll.2P n. hence 

.T . rt+h ^ 

l2< dt dx ds{Ae' -u' -Ve' --9%dwu'){s,x)-{w' *ja){t,x) 
Jo JR3 Jt 2 

(5.11) 

Then integrating by parts and by using (12. 6|) . with s = 0, p = oo and q = 2, 
we deduce 

/ ft+h 1 \ 1/2 

l2<h\\Ve^\\l2 +Ca~^/'^T^/'^\\e^\\LooL2 ih ||u"-Vr--(divn^)r lipids) 

< h\\vem,+Ca-^/^T'/^h\\n^LrLi\\u'-ve'- lidivu')e^LUi- 

(5.12) 

Summing up /i, I2 and by taking into account ()4.5|] . (14. 6p . (14. 7p . (|4.8p . we 
have 

||r(i + h)- ^"(i)||i2([o,T]xM3) < C(arV2 ^ /jc,-3/22^1/2 + /j)^ (5.13) 

by choosing a = h"^'^, we end up with (jS.Sp . D 
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As for 6^ we prove the equicontinuity in time for the sequence Pu^. 

Lemma 5.5. Let us consider the solution {u^,6^,p^) of the Cauchy problem 
for the system (jl.2p . Assume that the hypotheses (ID) hold. Then for all 
h G (0, 1), we have 

\\Pu'{t + h)- Pn^(t)||L2([o,T]xR3) < Ct/i^/^ (5.14) 

Proof. Let us set z^ = u^{t + h) — u^{t), we have 

\\Pu'{t + h)- Pn^(t)||i2([o,T]xM3) = 11^ dtdx{Pz') ■ {Pz' - Pz'*j^) 

+ 11 dtdxiPz') ■ {Pz' * j«) = h+ h. 

Jo JK3 

(5.15) 
By using ()2.5p we can estimate /i in the following way 

Ii < ar^/2||u^||Loo^2||VM=||r2 . (5.16) 

i- r^ II lli^t n^ II "'^t,x ^ ' 

Let us reformulate Pz^ in integral form by using the equation (II. 2P ]^. hence 



T r rt+h -j^ 



dt dxl ds{Au^ - {u^ ■V)u^ - -u%divu^){s,x) ■ {Pz" * ja)it,x) 
jr3 Jf 2 

(5.17) 



h< 



Then integrating by parts and by using ()2.6p . with s = 0, p = oo and q = 2, 
we deduce 

h < h\\Vu^\\% + Ca-^/^T^/^h\\u'\\LooL2 \\ {u' ■V)u' - -{dwu')u' 11^2^1 . 

t,x t X 2 t X 

(5.18) 
Summing up /i, I2 and by taking into account ()4.5p . (I4.6p . (14. 7p . ()4.8|) . we 
have 

||Pu^(t + /i) - i^u'(t)|li2([o,T]xR3) < C(ari/2 + /ia-3/2ri/2 + /i), (5.19) 
by choosing a = h"^'^, we end up with (|5.14p . D 

Corollary 5.6. Let us consider the solution {u'^,9^,p^) of the Cauchy prob- 
lem for the system (|1.2p . Assume that the hypotheses (ID) hold. Then as 

Pu' — > Pu, strongly in L'^{0, T; lI^{R^)). (5.20) 

r — > e, strongly in L'^{0, T; lI^{R^)). (5.21) 

Vr ^ Ve, weakly in ^^((0, T) x M^). (5.22) 

Proof. By using the Lemma 15. 5^ the Lemma 15.41 and the Theorem 15.31 and 
the Proposition 15.21 we get ()5.20p , ()5.2ip , while ()5.22p is a consequence of 

3D. n 
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5.2. Proof of the Theorem [331 



(i) It follows from the estimate (I4.6p . 
(ii) It is a consequence of the Proposition 15.21 

(iii) By taking into account the decomposition u^ = Pu^ + Qu^ , by the 
Corollary 15.61 and the Proposition 15.21 we have that 

Pu" — > u strongly in L'^{[0, T]; L^„^(M^)). 

(iv) Follows from the estimate (|4.6p and by using the Corollary 15.61 
(v) Let us apply the Leray projector Q to the equation ()2.1P ]^. then it 
follows 

V/ = AQu^ - Q ( div(ti^ » n^) + -u^ div Qn^ J . (5.23) 

Now by choosing a test function ip £ HlW^'"^^^ n C^Wx'^ n LfF^ 
and by taking into account ()4.5p . ()5.3p . (I5.20p . we get, as e | 0, 

(u^divQu%Q^) < ||Q'u"||L2^4||Vu^||i2i2||Q(/?||i-L4 

+ WQu'h^^Lih'hrLiW^Q^WL^Lt ^ o, (5.24) 



(div(n^ ^u^),Qip) = (div(Pu^ Pu^), Qip) + (div(Qn^ Qn^), Qv?) 

+ (div(P'u^ Qu^), Q(/7) + (div(Qu^ ® Qu^), Qv?) 

^ (div(Pu^Pn),Q(/7) = (Qdiv((Pn • V)Pn),(/7). 

(5.25) 

So as e I one has, 

(V/, if) — > (VA-^ div((u • V)u), 99). (5.26) 

(vi) Before passing into the limit in the system ()1.2p the convergence of 
the nonlinear term u'^-'VO'^ deserves a little discussion. By using again 
for u^ the associated Hodge decomposition, namely u^ = Pu^ + Qu^ 
and by taking into account ()5.3I1 , (I5.20p , (I5.2ip , (I5.22p one can prove 
that, for any p G [4, 6), 

Qu'-Ve' — >0 strongly in Li([0,r];L^(M3)), (5.27) 

Pu" • Vr — >u-V9 strongly in L^{[0,T] x M^). (5.28) 



Now we can pass into the limit inside the system (jl.2p and we get 
u and 9 satisfy the following equation in P'([0,T] x 



P{dtu - An + (n • V)u) = 0, (5.29) 

dt9-Ae + u-V9 = 0. (5.30) 

Finally we prove the energy inequality. By using the weak lower 
semicontinuity of the weak limits, the hypotheses (ID) and denoting 
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by X the weak- limit of \fE]f , we have 

^(|x(x,t)P + |^(x,t)|2 + |e(x,t)|2)dx 

JO JR3 
<hminf(- / (e|/(x,t)p + |n=(x,t)p + |r(x,t)p)dx 



e- 



^0 V2 



+ liminf( /" I (|Vn^(x,s)p + |vr(x,s)p)dxds 
= hminf / \ {e\v%{x)\^ + \ul{x)f + |0§(x)|2) dx 

^(|no(x)|2 + |0o(x)|')dx. (5.31) 
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